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Abstract 

The residual interaction for a meson-meson system is computed utilizing the cu- 
mulant, or cluster, expansion of the momentum-space time correlation matrix. The 
cumulant expansion serves to define asymptotic, or free, meson-meson operators. The 
definition of an effective interaction is then based on a comparison of the full (inter- 
acting) and the free (noninteracting) time correlation matrices. The proposed method, 
which may straightforwardly be transcribed to other hadron-hadron systems, here is 
applied to a simple 2+1 dimensional U(l) lattice gauge model tuned such that it is con- 
fining. Fermions are treated in the staggered scheme. The effective interaction exhibits 
a repulsive core and attraction at intermediate relative distances. These findings are 
consistent with an earlier study of the same model utilizing Liischer's method where 
scattering phase shifts are obtained directly. 
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1 Introduction 



Since the introduction of the 7r-meson by Yukawa [1] boson-exchange models have 
served to describe the phenomenology of the strong force, in particular the nucleon- 
nucleon interaction [2]. However, an understanding of the nature of the strong interac- 
tion between hadrons from first principles, namely quantum chromodynamics (QCD), 
still remains a central, but nevertheless elusive, problem in nuclear physics to this 
date. From the nuclear-physics perspective the existence of a confining phase in QCD 
is a defining feature since it is responsible for allowing only colour-neutral particles to 
propagate as asymptotic states. The mutual interaction of those particles is a central 
subject of nuclear physics. 

The physics of hadronic interaction probably is very complex. Apart from confine- 
ment, some prominent aspects to consider are the internal structure of the hadrons as 
determined by the quark-gluon- vacuum dynamics of the field theory, further, quark and 
gluon exchange, the Pauli principle, sea-quark loops, chiral symmetry, features of the 
SU(3) gauge group, quark flavors and masses, and more. Very likely, different mecha- 
nisms contributing to an effective residual interaction will do so in different dynamical 
situations. For example, sea-quark loops may dominate the long-range attractive part 
of the interaction, as those are related to virtual 7r-meson propagation in the language of 
boson-exchange models. At small relative distances quark exchange and gluon degrees 
of freedom may assume the dominant role. 

Clearly, a study of hadronic effective interaction based on first principles needs to 
implement the nonperturbative aspects of the quantum gauge field theory. This is cur- 
rently best done using lattice gauge field theory and numerical simulation. The number 
of attempts to understand the nucleon-nucleon interaction from lattice QCD are few. 
The work of Rabitsch et. al. [3] employs the finite-temperature formalism, no repul- 
sive core is seen. Further, despite huge computational obstacles, recently progress has 
been made towards the qualitative understanding of 7T-N and N-N scattering lengths in 
QCD with heavy quarks and gluon degrees of freedom [4, 5]. Since the latter simulation 
already pushes contemporary numerical resources to their limits, the question arises 
which avenues can be pursued to study the physics of effective hadronic interaction 
beyond the zero-momentum limit. 

The philosophy of the present work is that one may use other, simpler than QCD, 
field models to capture the basic physics of the effective interaction. The minimal re- 
quirements towards such a lattice gauge model is that it is confining, at least for certain 
choices of the coupling, and that it contains fermions. The hope is that essential dy- 
namical features pertaining to the effective hadronic interaction will be qualitatively 
similar in QCD. Another important aspect is the development of theoretical concepts 
and computational techniques which may then be applied to the QCD case in subse- 
quent simulations. 

In a previous work [6] meson-meson scattering phase shifts were obtained in a 2+1 
dimensional lattice model with a U(l) gauge group and staggered fermions (QED 2+1 ). 
This model evidently oversimplifies structural details of QCD 3+1 . For instance, due to 
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the availability of only one colour baryons are absent, the hadrons in the theory are 
mesons. On the other hand QED 2+1 posesses a minimal set of features which should 
be present in any model aiming at effective hadronic interactions: 

• Confinement, at least for some range of the gauge coupling, which allows the 
propagation of asymptotic colour-singlet particles. 

• Fermions, so that a set of hadron-hadron operators can be constructed. 

• Space dimension d > 2, so that various partial waves, or spin channels, can be 
studied. 

Indeed, the results of the scattering simulation [6] were consistent with the gross fea- 
tures of the N-N interaction (short-range repulsion, intermediate-range attraction). On 
the other hand, Liischer's method [7], which was utilized to obtain the phase shifts, 
possesses inherent ambiguities which make it very difficult to interpret the simulation 
data. Additional, a-priori, assumptions about the effective interaction are almost a ne- 
cessity. Furthermore, phase shifts are only obtained for a discrete, finite, and typically 
very small, set of momenta. Interpolation between those is only possible by repeating 
the simulation on a variety of different-sized lattices. 

It is the purpose of this article to develop a viable alternative to Liischer's method. 
We will subsequently describe how an effective hadron-hadron interaction can be de- 
fined within the framework of a lattice field theory, and how it can be extracted from a 
numerical simulation. Although the presentation proceeds within the QED 2+1 model, 
generalization to other lattice models is straightforward. It is understood (though not 
attempted here) that the resulting effective interaction may eventually be employed to 
calculate scattering phase shifts and other two-hadron properties. 

Finally, in judging the numerical quality and shortcomings of the present simulation, 
the reader should be aware that lattice nuclear physics, which we may loosely define as 
the "chemistry of colour-neutral objects on the lattice" is subject to a change of scale: 
The residual interaction energies are only 10~ 2 -10~ 3 times those of a typical hadron 
mass. Nevertheless, the proposed method, which involves taking ratios of simulation 
data, is capable of dealing with the transition from the GeV to the MeV scale, although 
at reduced numerical quality. 

2 Lattice Model 

As in ref. [6] we use an L 2 X T = 24 2 X 32 lattice with the compact Wilson action for 
a U(l) gauge group 

S w =(3j2(l-ReU D ) (1) 

□ 

where Un = e @n and ©□ is the oriented plaquette angle. For the fermions the Kogut- 
Susskind action 

Sf=J2 Xf{x)G-l[U} Xf {y) (2) 
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is employed [8, 9]. The Grassmann fields % and % have been given an external flavour 
index / = u,d and 

G x?y = ^ E Vni x ) [Un( x )6x+vL, y ~ Ul(y)6 x , y+I1 ] + m F 8 XtV (3) 

is the staggered fermion matrix with link variables U^(x) and fermion mass rap. Units 
are such that the lattice constant is a = 1. At (3 = 1.5 the lattice model is well in the 
confined phase, see [6] for details. 

3 Meson Fields 

One-meson fields are constructed from boosted local operators 

fat) = L- 2 £ j**Xd& t) Xu (x, t) (4) 

X 

where the sum extends over the L 2 sites of the spacial lattice. Thus describes a 
meson (resembling the ir + ) with momentum 

p=^-(h,k 2 ) where k 1>2 = -(^ - 1) . . . ^ , even L . (5) 

Two-meson fields with total momentum P = then are 

^{t) = <j ) _ ? {t)<j )+? {t). (6) 

Correlations of these operators contain the information about the dynamics of the 
meson-meson system and, ultimately, the effective residual interaction. 

4 Time Correlation Matrices 

The 2-point correlator, describing the propagation of one meson on the lattice, is 

c$(Mo) = <4w^°)> - <4w><^°)> • (7) 

It can be worked out in terms of contractions between the Grassmann fields, say 

n n J n 

■ ■ ■ *f ( x ) ■ ■ ■ Xf> ( x ) • • • = • • • Sff Gxx' ■■■ (8) 

where n indicates the partners of the contraction, and G is the inverse fermion matrix, 
see (3). We obtain 



c#(Mo) = i- 4 EE e " ip - x+l9 - y (l^ol 2 ) • (9) 

x y 
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Assuming translational invariance in the form 

(\Gg t ,gt 1 2 } = (\Gg + a t ,y+at 1 2 } (10) 
one of the site sums in (9) can be worked out, rendering the 2-point correlator diagonal 

, *o) = S M L- 2 £ e-* g (\G xX z oto | 2 } . (11) 

X 

The point x is arbitrary, with the correlator being independent of x . We will later 
describe how is computed and explain its importance for extracting the effective 
interaction. 

The 4-point correlator describes the propagation of two interacting mesons on the 
lattice 

C$(Mo) = <<#)^°)> " W> • (12) 

Here p and q are the relative momenta in the meson-meson system. In terms of fermion 
propagators this correlator reads 



c«(Mo) = l- 8 J2 E E E e^-^+w-*: 
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(13) 

which we here list for the sake of completeness. More important for the present purpose 
is a careful analysis of with regard to the meson-meson effective interaction which 
is somehow contained in it. Towards this end a diagramatic classification of the various 
terms contributing to proves useful. Such a classification arises naturally from 

working out the contractions between the quark fields in (12), (6), (4). Let us write 

= 1 1 + X - M - M < 15 > 

for the four terms as they occur in (13). Each diagram line corresponds to a fermion 
propagator. For example, using the notation introduced in (8), we have 

43 21 12 34 43 34 21 12 

C {4A) = <^V_ f * f > = <^V_ f > (16) 

21 43 12 34 43 34 21 12 

C (4B) = <^V_ f * f > = <*W^-f> (17) 

where the n = 1 ... 4 identify the partners % and % of a contraction, see (4) and (8). 
Clearly, in diagrams A and B all pairs of contractions preserve the integrity of the 



one-meson fields cj) and ft in the sense that no breakup via quark exchange occurs. 
At the most, diagram B describes the exchange of (composite) mesons as a whole. As 
we proceed, diagrams A and B will have to be treated on the same footing since the 
desired effective interaction needs to respect Bose symmetry on the meson level. 

Diagrams C and D on the other hand both contain quark exchange between the 
mesons 

23 41 12 34 

C (4C) = <^V_ f * f > (18) 

41 23 12 34 

C {iD) = <4^-^> (19) 

41 i 4 

for example ft ... (fx/) etc. Thus these diagrams must be considered solely as sources 
of the effective interaction. 

It should be emphasized that the propagation of free, noninteracting, mesons is not 
quite described by C^ 4A ^ + C^ 4B \ The reason is that gluonic correlations are introduced 
between the mesons by taking the gauge configuration average (. . .) over the product 
of all four operators. (This is of course also true for C^ 4C ^ and C^ 4D \) Nevertheless, it is 
clear that the uncorrelated propagation of two mesons on the lattice is fully contained in 
(ji iA ) _|_ C^ 4B \ It is essential to isolate the uncorrelated component because it identifies 
the channel, in the sense of scattering theory, relative to which an effective interaction 
can be defined. 



5 Uncorrelated Mesons 

In order to identify the free, noninteracting, part of the meson-meson time correlation 
matrix let us, for example, examine C^ 4A \ Expressing the contractions in (16) in terms 
of quark propagators, with the notation introduced in (8), we have 

C( 4 ^) ~ (G*G G*G ) (20) 

where the ~ indicates the Fourier sums, etc., which carry over from (4). For full math- 
ematical detail see (13). 

The gauge configuration average in (20) may be analysed by means of the cumu- 
lant (or cluster) expansion theorem, see for example section 12.3 of ref. [10]. Taking 
advantage of (G) =0 we have 

(G*G G*G ) = (G*G )(G*G } + (G*G )(G*G } + (G*G*)(G G ) + ((G*G G*G )) . 

(21) 

The last term defines the cumulant (which generalizes the standard deviation). The 
first three terms on the right-hand side of (21) are illustrated in Fig. 1. The dashed 
lines indicate that the propagators are correlated through gluons. Evidently only the 
first one of the three separable terms in (21) represents free, uncorrelated, mesons. All 
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Figure 1: Illustration of the cumulant expansion (21) of the 4-point correlator. 



other terms in (21) are sources of residual effective interaction between the mesons. We 
therefore define 

43 34 21 12 43 2 1 

C [4A) = <^ f WV_ f > ~ (G*G )(G*G ) . (22) 
A similar analysis of C^ 4B ^ leads to 



43 34 21 12 



^ (4B) = (^)(4^)- ( 23 ) 
The sum of those then is the free meson-meson time correlation matrix 

t ) = C% A \t, t ) + C% B \t, to) (24) 

which describes two noninteracting mesons on the lattice. We emphasize that is 
an additive part of the full 4-point correlator. 

C« = <7« + C} 4) (25) 



The remainder C\ contains all sources of the residual interaction, be it from glu- 
onic correlations, or quark exchange (or quark-antiquark loops, if the simluation is 
unquenched) or any other complicated effects. 

It also should be emphasized that the free correlator describes (isolated) mesons 
the masses of which, as well as other internal structure features, consistently arise from 
the dynamics determined by the lattice field model and its numerical implementation. 
From the point of view of numerical simulation it is crucial to capture the small (10~ 2 - 
10~ 3 ) effects of the residual interaction relative to the computed free correlator. 

In the light of the above, the free meson-meson correlator should be obtained 
numerically from the quark propagator. A glance at (22), (23) and (7) suggests that 
may be expressed in terms of the 2-point correlator C^ 2 \ We note in passing that 
the separable term in (7) is identically zero because of the u, d flavor assignment to the 
Grassmann fields of the one-meson operators (4). Thus, using (7) and (22), (23) gives 

C£ = C%CW_,+ C%C^. (26) 

Also, recall that is diagonal in the momentum indices. Writing 

Ci 2 )(t,t ) = 6 e? c P it,t ) , (27) 

where c^(t,t ) is given by comparing (27) and (11), and has the obvious property 

c_ (t,t o ) = c}(t,t ) , (28) 

we obtain 

C$(t,t ) = \c P it,t )\ 2 . (29) 

In this form is easily computed, for example using (9) to obtain c^(t,t ) = 

Cl^(t, t ). This can be done in a manner consistent with the computation of Cl^(t, to), 
for example using (13). Another property, which will become important later, is that 
the eigenvalues of are proportional to \cp(t,t )\ 2 and thus are nonnegative. Finally, 
it is interesting to note that the property 

Oi 4 l = C (4 h = 0i 4) , = C {i l , (30) 

p,q -p,q p,-q ~P,~q V J 

which is evident from (26) and also holds for C^ 4 \ reflects Bose symmetry with respect 
to the (composite) mesons. Permutation of the mesons, one with momentum -\-p the 
other with momentum — p, results in the substitution p — > —p. Clearly the symmetry 
(30) stems from the fact that was defined through the sum of diagrams C^ 4A ^ and 
C( 4B \ see (24). 



6 Effective Interaction 

In principle, the deviation of the full time-correlation matrix from the free time- 
correlation matrix contains all the information about the effective meson-meson 
interaction. 
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A rigorous definition of the latter emerges from exploring the properties of a time 
correlation matrix for an elementary interacting meson field subject to canonical quan- 
tization. This is discussed in appendix A. In first-order perturbation theory an ex- 
plicit formula for the (small) interaction can be derived. It involves the zeroth-order 
correlator, which represents the noninteracting system, and the full correlator, which 
represents the interacting system. By way of analogy we adopt this formula to guide 
us towards defining an effective interaction in the framework of a lattice simulation. 

Thus, in view of (68), let us define the effective meson-meson time-correlation matrix 



C«(Mo) = <7M(Mo)" 1/2 cM(Mo)<? (4) (Mo 



-1/2 



(31) 



As explained in appendix A multiplying by the two inverse square roots of C^ 4 \ 
in a symmetric way, neutralizes the effects of compositeness of the individual mesons, 
as well as the trivial aspects of free propagation. The effective correlator deviates 
from the unit matrix (constant for all t) only if a relative residual interaction between 
the mesons is present. The latter then is defined as 



This is obviously consistent with 

C (4) (Mc 



<9C( 4 )(Mo) 



dt 



(32) 



t=to 



-H J (t-t ) 



(33) 



In appendix A we verify the validity of (the elementary meson field version of) formulas 
(32) and (33) up to first-order peturbation theory. Equation (33), in conjunction with 
(31), may simply be adopted as the definition of a residual effective interaction. On the 
other hand, perturbation theory should be perfectly valid since the effective interaction, 
as it comes out of a lattice simulation, is about 2-3 orders of magnitude smaller than 
a typical hadron mass. 



7 Numerical Procedures 

Computation of the 4-point correlation matrix using the explicit form (13) is nu- 
merically not feasible. Not only would it be necessary to compute the entire fermion 
propagator matrix Gst,yt , but also, the four-fold sum over the spacial lattice contains 
L 8 = 10 11 terms. This computational problem is solved by using a random source 
technique adapted to the Fourier-type sums which occur in (13). We have used ran- 
dom sources exactly as described in [6], including the choices for all parameters. All 
correlation matrices were computed in momentum space with a truncation given by 

p=^-(k 1} k 2 ) with 1*1,2! < 2. (34) 
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The 2-point correlator (which is not needed in Liischer's method) is obtained from 

Cm = <E E ^ t , to (0; i2)W?i R)) ■ (35) 



P>9 

<i?> 



Here Hg tito (p; R) are propagator columns corresponding to Fourier-modified Gaussian 
random sources exactly as defined in equation (26) of ref. [6], and S<.r> denotes the 
average over Nr = 32 random sources R. Equation (35) is particularly useful for testing 

(2} 

the diagonality as claimed in (11). We have verified C~~(t,to) oc 8$$ for all momenta 
allowed by (34), within statistical errors. 

Since both the full and the free correlators and respectively, commute with 
all elements of the cubic lattice symmetry group 0(d, Z) in d = 2 (space) dimensions 
one may compute their reduced matrices, say C^ 4 ' r ^ and C^ 4 ' r \ within each irreducible 
representation T of 0(d, Z). This has been discussed in detail in ref. [6]. We here just 
mention that there are five irreducible representations T of 0(d, Z), four of which are 
one-dimensional T = A\ , A 2 , B\ , B 2 , and one is two-dimensional T = E. The latter is 
excluded on grounds of the meson-meson system having positive parity. The basis states 
within each of the former representations spaces T can be labelled by the magnitude 
p = \p | of the momenta. This is unique for the truncated set (34) of momenta. As is 
obvious from (29) the reduced free correlator then has the form 



3£ r) (Mo) = S pq 4 r >(Mo) 2 (36) 

for T = A\ , A 2 , B\ , B 2 , where c^(t,t ) are linear combinations of 2-point correlator 

elements CiJ, with p' = O g p, g (E 0(d, Z). Those can be easily constructed with the 
transformations given explicitly in appendix C of ref. [6] . 

The time behaviour is dominated by an exponential function at large t. Since we 
employ periodic boundary conditions this means 

cf>(Mo)| 2 = 4 r) cosh(wf>(*-g) for t^t c (37) 

with t c = T/2 + 1 = 17, recall T = 32. The strength factors A^ and the free two- 
meson energies Wljf^ are listed in Tab. 1. They were obtained from least-square two- 
parameter fits to the data using (37). All errors throughout the present work stem 
from a jackknife procedure [11] with Nu = 64 omissions. Some of the momenta p 
belong to more than one sector T, as indicated in Tab. 1. This degeneracy is realized 
analytically. The strength factors become very small for large relative momenta p. This 
reflects unfavourably on the errors in the last line of Tab. 1. Eventually, data which 
derive from the largest momentum index, k 2 = 8, will be dropped from the subsequent 
analysis. 

Given (36) it is trivial to construct the inverse square roots of the 2-point correlation 
matrices which occur in the definition (31) of the effective correlator. The reduced 
matrix elements of the latter thus are 

nWft i \ - ^/^(Mo) , . 
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Table 1: Strength factors and energies for two noninteracting mesons obtained from 
the free corrrelator |c£ r )(t, t )\ 2 according to (37). The momenta are referenced using 
(34). 

The effective correlator matrix elements are obtained as a ratio and as such are directly 
sensitive to the sought-after subtle effects of the residual interaction. 

In order to obtain the time derivative of the effective correlator matrix C^ 4 ' r \t, t ), 
see (32), we have first numerically diagonalized the full correlator, say 

AT(r) 

Cif\tM = E vPWW^vPiq) (39) 

n=l 

with vP(p) being the components of an orthogonal system of real unit vectors 

vP within each of the four allowed sectors T. At large t, again meaning t = t c , the 
eigenvalues behave like 

\P(t,t ) = A£) cosh (wP(t-t c )) for t = t c (40) 

whereas the eigenvectors are independent of t. Time-independence of the eigenvec- 
tors has been verified numerically for five time slices, t = 15 . . . 19, around t c = 17, by 
computing the statistical errors on the components of the . 

The values for A^ and WP in (40) were obtained from two-parameter fits to the 
data using five time slices 2 = 15... 19. Since this has also been done in the previous 
work [6], where the same gauge field configuration were used, the results are identical 
with those in table 2 of ref. [6]. In the sectors A\ and B 2 the smallest eigenvalues 
\([\t,to) for n = 6 and n = 3 respectively (which belong to the largest energies) 
are contaminated by round-off errors beyond control. As mentioned earlier, we face a 
similar situation in the analysis of the free correlator, as exemplified by the last line of 
Tab. 1. Thus, our final choice for the sizes of the correlation matrices C^ 4 ' T \ C^ 4 ' r ^ and 
C( 4; H was ]y( r ) = 5 ; 1 ; 3 ; 2 for the sectors T = A\, A 2 , B\, B 2 respectively. 

In order to reconcile the cosh time behavior of the correlators with the pure exp 
behavior assumed in appendix A, and (33), we may choose t — t c large and negative 
in (40) and (37), and then drop the negligible ones of the exp terms in the hyperbolic 
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cosines. In this way, from (36)-(40) we obtain 



AT( r > .(r) 
4f } (Mc) = E V P(P) -r===z= ex P 

n=l 



' , m py(r) + p7(r)\ 



*< r) (?) 



The time derivative in (32), which now has to be taken at t = t c , then leads to 



(41) 



These are the desired matrix elements of the effective interaction in the basis \(T;p 2 )). 
From a numerical point of view the computation of those matrix elements is the most 
critical step of the simulation because it is in (42) where subtle effects of the effective 
interaction reveal themselves. 

To illustrate the numerical situation we have listed in Tab. 2 the components v^(p) 
of the eigenvectors, some strength ratios 



(P. 9) 



and energy shifts 



wp + wp 



(43) 



(44) 



both for p = q, and for the sector r = A\. As might be expected, given the smallness 
of the residual interaction, the effective correlator is diagonal to a very high degree. 
This is reflected in the observation that v^\p) % 1 for only one component of 
with all others being quite small. In all cases, however, the largest component of 
is multiplied, in (42), with an energy shift w£\p,p) of zero, within error limits, see 
Tab. 2. Also, the strength ratios a^(p,p) increase as p becomes larger. This is true for 
all n, however, their overall size is largest for small n. It appears that the contributions 
to the momentum-space matrix elements are mostly driven by the behaviour of the 
strength ratios, as these vary most dramatically, with their "direction" determined by 
the signs of the energy shifts. Off-diagonal contributions carry factors v^(p)v^(q), 
with p ^ q, and are thus somewhat supressed. Those are, however, not negligible. 
In summary, the numerical situation is best characterized by saying that the matrix 
elements {p 2 \H { p\q 2 ) emerge in a subtle way from the various factors and terms in 
(42). The other sectors r = A\, B\, B 2 behave similarly, although the magnitude of the 
factors is significantly smaller. 

In Tab. 3 we show the momentum-space matrix elements (p 2 \'Hj \q 2 ) for the sector 
r = A\. Diagonal elements are numerically zero. The matrix elements are largest in 
the off-diagonal corners of the matrix. 
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—0 01 (61 




5 


-0.02(7) 


3.23(5) 


-0.05(6) 


5 





-0.013(1) 


0.0010(2) 


0.38(9) 




1 


-0.026(3) 


0.0049(9) 


0.30(9) 




2 


-0.026(9) 


0.022(4) 


0.21(9) 




4 


0.02(7) 


0.29(5) 


0.03(9) 




5 


0.999(1) 


0.9(2) 


0.00(9) 



Table 2: Selected factors contributing to the matrix elements of the effective interaction 
according to (42)-(44). The momenta are referenced like in Tab. 1. 








1 


2 


4 


5 





-0.01(4) 


0.004(3) 


0.009(4) 


0.029(8) 


0.08(3) 


1 




-0.02(4) 


0.009(5) 


0.021(7) 


0.05(2) 


2 






-0.03(4) 


0.004(4) 


-0.000(6) 


4 








-0.03(6) 


-0.05(1) 


5 










-0.1(1) 



Table 3: Matrix elements (p 2 \7ij \q 2 ) for the sector T = A\. Row and column indices 
are labelled by k 2 , as in Tab. 2. 
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It should be noted that numerical alternatives to the preceding analysis have been 
explored. For example, fits to c^(t, t ) with a cosh function like (37) lead to one-meson 
energies which should be close to Wljf^ /2. However, given that 

COsh[W(t - t c )] = COsh 2 [y(t - t c )] + Sinh 2 [y(t - t c )] , 

this is the case only to some degree of approximation because of contamination from 
the sinh term. Since the definition of the effective interaction is based on a comparison 
of two meson-meson correlators, and C^ 4 \ it is more consistent to fit to data from 
C^ 4 \ as put forward in (36)— (37). 

Another alternative is to compute the effective correlator £( 4;r ) directly from the 
ratios (38), then diagonalize it, and make cosh fits to the eigenvalues with, say 

«i r ) cosh[ e (P( t - Q] . 

This will yield in the diagonal representation, with eigenvalues oc ct^ \ \ . This 
approach appears very appealing because the ratio data (38) are taken advantage of 
directly. Unfortunately, the errors on the matrix elements C^ r \t,t ) tend to be am- 
plified by the diagonalization routine to a larger degree than those of the interacting 
correlator matrix elements C^ r \t,t ). This effect originates with the ratios (38) for 
the larger momenta p and q where the numerical values of the numerator and denomi- 
nator become quite small. The effective interaction obtained in this way appears to be 
consistent with the results shown below, however, the jackknife analysis of the errors 
proves difficult because for some of the samples an acceptable (small % 2 ) cosh-type fit 
could not be obtained. 

The results presented in the next section are based on employing (42). 

8 Coordinate- Space Effective Interaction 

Given the matrix elements (42) it is straightforward to construct the coordinate-space 
representation of Tii 

{t\Hi\?') = E( f lM n IO (45) 
r 

= E E (r\(T;p 2 ))(p 2 \nP\q 2 )((T; q 2 )\r') . (46) 

r p 2 , q 2 

The transformation matrices (f\(T;p 2 )) between the basis | (T; p 2 )) of the sector T 
and lattice coordinate points \r ) are constructed by applying appropriate projection 
operators to (r \p) = L~ 1 e w ' r . We refer the reader to appendix C of ref. [6]. 

Of the four representations T = A\ , A 2 , B\ , B 2 which contribute to (r \Hi\r') the 
A\ sector numerically dominates the T sum in (45). The Ai-sector contribution to the 
diagonal elements (f\Hi\r) versus r = \r \ is shown in Fig. 2. One source of fluctuations 
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Figure 2: Contribution of the A-y sector to the diagonal matrix elements (r \hv I \r) of 
the effective interaction versus r = \r |, according to (45). 
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Figure 3: As Fig. 2 for the A 2 sector. 
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of the data points is the lack of continuous rotational symmetry on the lattice. For 
example, matrix elements (r\l~6^\r), in V = A\, are different for f = (5,0) and 
r = (4, 3), whereas r = 5 in both cases. This effect is more pronounced for larger values 
of r, and significantly more in the other sectors Y . For illustration Fig. 3 shows the A 2 
sector contributions to (45). 

Since the spacial extent of the lattice is L = 24, the maximal relative distance 
between two mesons, in the presence of periodic boundary conditions, is r = 12 X \/2 = 
17 with r along a diagonal. However, the data tend to repeat themselves already starting 
at r = 12 which corresponds to realigned along one of the axes of the lattice. Multiple 
interactions "around the world" between mesons in different copies of the lattice should 
be expected for 12 < r < 17. Interpreting the data in the context of an effective meson- 
meson interaction is therefore meaningful only in the region r < 12. 

Statistical errors in Fig. 2 and Fig. 3 stem from a jackknife analysis with Nu = 64 
omissions. Those matrix elements {p 2 \H^\q 2 ) which were less in magnitude than their 
statistical errors have been omitted in the transformation (46). 

We have tested the stability of the effective interaction by truncating the matrix 
in Tab. 3 down to size 4x4. The coordinate-space effective interaction, which is dis- 
cussed below, did change by about one cr, but was still within error limits. It should 
not be surprising that the qualitative appearance is not changed too much by the trun- 
cation since the large-momentum components are only responsible for high-resolution 
features. The test also shows that the small-momentum approximation to the effective 
interaction has some degree of justification. 

The A 2 -sector contribution to (r \Hi\r') exhibits a very sensitive dependence on 
the direction of r, as is evident from Fig. 3. This is also true for the sectors B\ and 
B 2 (not shown), but to a much lesser extent for the A\ sector, see Fig. 2. A sensitive 
angular dependence is indicative of the presence of nonzero angular momenta. We have 
used a projection technique, described in appendix B, to interpolate our lattice data for 
(f\Hi\r') with an 0(2) rotationally invariant operator. Its partial- wave decomposition 
is given by the im-sum in 

(r \Hj\f ') = £ Y lm ^)H?(ry)Y; m ^') + . . . (47) 

with r = \r |, r' = |r'|, and <£>,<^' being the polar angles of r and r' respectively. 
Here the Yi m (tp) oc e %mt f are (circular) harmonics, as defined in (71) of appendix B. 
The sum extends over even angular momenta I = 0,2,4. . . only. In appendix B we 
elaborate on how the partial-wave I interpolants 7i\ (r, r') to the simulation data of 
the effective interaction are constructed (see (76) for a quick orientation). The ellipses 
in (47) indicates (small) terms which remain after angular momentum projection, and 
which in principle should vanish in the continuous limit (a — > 0, L — > oo) as the 
discrete 0(2, Z) symmetry approaches continuous 0(2) symmetry. A good indication 
of the smallness of those terms, in partial wave I = 0, is provided by comparing Fig. 4, 
which shows Ti^\r,r) for I = 0, with the "raw" data of Fig. 2. It appears that "cubic" 
fluctuations are simply smoothed out. In partial waves 1 = 2 and I = 4 the "raw" data 
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Figure 4: Diagonal elements T~6p(r,r) of the angular-momentum-projected effective 
interaction, for partial wave I = 0. 
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(for example see Fig. 3) are weighted by factors with alternating signs upon projection. 
This results in very small 1 = 2 and I = 4 wave components for the effective interaction. 

As can be seen from (76) and Tab. 4 of appendix B, the partial-wave projected 
effective interaction T-6p{r,r') receives the following contributions: 

I = 0: A-l only 
1 = 2: B t and B 2 
£ = A: A-y and A 2 . 

The simulation data clearly show that the partial wave I = dominates over 1 = 2 
and I = 4. 

The effective meson-meson interaction Tii obtained from the present lattice simu- 
lation is non-local. This is an inevitable consequence of the momentum-space approach 
used in the present work, see (46). In Figs. 5,6,7 we show 3D plots of T-6j\r,r') for 
I = 0,2,4 respectively. 

The salient features of Ti.^~°\r , r') are repulsion at short distances and a region 
of attraction at intermediate relative separations. These properties are consistent with 
earlier results from a scattering phase shift simulation for the same system [6]. For 1 = 2 
the interaction is much smaller, see Fig. 6, with a repulsive soft core at intermediate 
r, r' and some attraction at large distances. As I increases, the structures are naturally 
pushed out to larger r,r' as the threshold behavior (~ r l ) becomes more dominant. 
This has apparently happened for I = 4 to an extent that makes it unfeasible to observe 
a possible attractive feature at large r, r' because, as mentioned above, r < 12 is the 
upper limit for a valid interpretation of Tii as a meson-meson interaction. 

It is important to be aware of the limitations of the present results. The small- 
momentum approximation, see (34), naturally limits the resolution of the method with 
respect to small structures in 7i\ (r, r'). The largest momentum used in the simulation, 
see Tab. 3, corresponds to a quarter wave length A/4 of about 4 lattice units a. In 
particular, the repulsive core, which is a prominent feature in Fig. 5, inevitably appears 
soft and broad as a result of the Fourier-type approximation. The latter also determines 
the "visual flavor" of the figures. 

It is of some interest to pursue an alternative approach based on a simulation 
aiming directly at coordinate-space matrix elements of Tii. In some sense the latter is a 
complementary (and incompatible) approximation which has drawbacks of its own, but 
nevertheless preliminary results [12] seem to indicate a much narrower (hard) repulsive 
core than Fig. 5 would suggest. A discussion of the latter approach will be forthcoming 

[13]- 

In order to examine the effective interaction in terms of, for example, gluon-exchange, 
quark-exchange, etc., the simulation data will have to be analysed in various other ways. 
Some of these aspects are currently being targeted by on-going work. However, the ori- 
gin of the repulsive core (in either of the two approximation schemes) clearly is the 
anticommuting (Grassmann) nature of the quark fields which imprint certain proper- 
ties [13] on the effective correlator, at least for the presently studied simple lattice field 
model and the corresponding class of meson operators. 
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Figure 7: Nonlocal g-wave interaction, I = 4. 
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9 Summary and Conclusion 



A method has been developed which allows to extract an effective hadron-hadron in- 
teraction from a lattice simulation of a gauge field theory with fermions. Although the 
practical application has, for now, been limited to a simple field model, a U(l) invariant 
confining theory in 2 + 1 dimensions with staggered fermions, the method is eventu- 
ally intended for studying hadron-hadron interaction in QCD 3+1 . The present proposal 
also is intended to provide an alternative to Liischer's approach [7] which suffers from 
inherent ambiguities related to the interpretation of scattering phase shifts. 

In the present application we were able to extract a residual effective meson-meson 
interaction Tii which features short-range (Pauli) repulsion and intermediate-range 
attraction, thus confirming the interpretation of earlier results [6] using Liischer's ap- 
proach. It should be emphasized that the effective interaction Tii is about one to two 
orders of magnitude smaller than a typical hadron mass. Extracting Mev effects from 
numbers on a Gev scale represents an extraordinary numerical challenge. Neverthe- 
less, it has been demonstrated that this program is feasible in a lattice simulation, 
although, naturally, the numerical quality that can be achieved is significantly less 
than in a standard hadron mass simulation, for example. 

In principle the proposed method straightforwardly adapts to applications with 
QCD 3+1 . Without doubt, problems of numerical magnitude will present themselves. 
Since it was possible to perform the present QED 2+1 simulation on a RISC Workstation, 
there is hope that the use of more powerful computing resources will make simulations 
of effective hadronic interaction feasible. 

It is a pleasure to thank H. Markum for interesting and useful discussions. We also 
wish to acknowledge the hospitality of CEBAF during the final phase of this work. 
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10 Appendix 



A First-order Perturbative Interaction 

We here elaborate on the definition of the effective meson-meson correlation matrix (31) 
and the effective meson-meson interaction (32) to the extent that those will be derived 
in first-order perturbation theory for the analogous case of an elementary interacting 
boson field. 

Let C = £o((j), d(j>) be the free lagrangian for an elementary boson field (j>{x) defined 
on the sites x = (x,t) of the lattice. It is understood that <j> is subject to the usual 
canonical quantization, say through commutators. Be £ = £ -\-£i such that Cj = Ci((j)) 
is a (small) interaction. In the usual way £ gives rise to a hamiltonian 

H = H + Hi (48) 

where H is the free part and Hi a perturbative interaction. With view on (4) and (6) 
define 

Mt) = L- 2 E e**#(d?, t) (49) 

X 

and 

= (50) 

In the correlation matrix 

c^%Q = (o\H(t)K<Q\o) (si) 

the separable term, see (12), has been dropped since it is zero for the flavoured quark 
fields, see section 4. The nondegenerate vacuum state |0) satisfies H\0) = Wo\0). We 
will assume that its energy is zero, Wo = 0. Thus the time dependence of the correlator 
(51) may be made explicit 

C$(Mc) = (0|$t(t c ) e -^-^^(t c )|0). (52) 
Switching to the interaction picture, define 

Hft) = e #o(t-t e ) fy i e -H (t-t c ) ( 53 ) 

and the (euclidean) time evolution operator 

U(t,t c ) = e Ho(t~tc) e -H(t-t c ) _ ( M ) 

The perturbative expansion of the latter 

oo (_i)N . t . t 
U{t,tc)= E / dtl '-- dt N T[H I {t 1 )...H I {t N )} (55) 
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then induces a perturbative expansion of the correlator 



Ci%Q = (0\^(t c )e- 6 ^U(t,t c )^(t c )\0) (56) 

oo 

= E^lMc). (57) 

N=0 

The zero-order term (7( 4 ' Ar=0 ) evidently describes noninteracting mesons. 
Order N = 

Let \nv) be a complete orthogonal set of eigenstates of H 

H \nu) = W^\nv) (58) 

where are free two-meson energies on the lattice, and v is a degeneracy index. 

Lattice effects set aside, those energies should be close to 2-^/m 2 + p 2 with m being the 
rest mass of one meson and p = \p | its momentum. Also, define the relative meson- 
meson momentum-space wave functions ip^)(p) through 

c^MJ(p) = <H^(*c)|0}* (59) 
where c^J are normalization factors. The order N = correlator then is 

^r =0) (*. *c) = E \^\ 2 e~ W ^-^Kv WW (<f ) • (60) 

With properly chosen normalization factors c^J we expect orthonormality and com- 
pleteness 

E ^'(p^fflCp) = <W>, (6i) 



V 



E^SCWS'C? ) = ( 62 ) 

For a free elementary boson field this is almost a trivial point since the ip^J merely are 
plane (lattice) waves. A glance at (60) shows that, technically, those could be obtained 
as (normalized) eigenvectors from diagonalizing the correlation matrix (j( 4 '> N = )^ where 
|c^| 2 are the eigenvalues, at t = t c . For the case considered here all eigenvalues of 

£(4;AT=0) wiU be nonzero . 

Considering a lattice model, however, where the role of <Ey(t c ) is assumed by a 
composite operator made from fermion fields, see (4) and (6), it can not be a priori 
excluded that an operator matrix element, of the type as it occurs in (59), is identically 
zero for all p. In this case the free correlator of the lattice model would have an 
eigenvalue zero (for all t). Likewise, if, for some reason, the set of hadron operators 
used to construct the correlation matrices on the lattice is linerarly dependent, one 
should expect to have an eigenvalue zero. 
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Order N = 1 



From (55)— (57) we obtain 

Cjr =1) (M«0 = -mi(t c )e- A ^ fdUH^^iQlO) . (63) 

J tc 

Upon inserting the complete set \nv) on both sides of Hj{ti) and using (53) the t\ 
integral over exponentials can be carried out explicitly. The result is 



a 



(4;AT=i; 



PI 



(t,t c 



( nv \ H i\m^) c^c^l exp 
sinh 



ti. 1 r r m 



2 



w i0) -w i0) 



(64) 



Without loss of generality the normalization constants c^J may be choosen real and 
positive, with the phase factors being absorbed into ip^)(p), as is evident from (59). 
Thus a glance at (60) shows that the two normalization factors and the exponential in 
(64) may be removed by multiplying the correlation matrix C( 4 > N=1 ) from both sides 
with the inverse square root of (j( 4 '> N = ). Hence the matrix elements of 



c^ N=i \t,t c ) = c^ N=o \t,t c y 1/2 c^ N=i \t,t c ) c^ N=o \t,t 



-1/2 



(65) 



in the basis ip^J{p) are products of (nv\Hi\m[i) and the expression inside {. . .} of (64). 
The t derivative of the latter is equal to one at t = t c . Thus we have 



dC, 



(4;AT=i; 



PI 



(t,tc 



dt 



EE^°J(^(^|^|m^t(f 



(66) 



t = tr 



nv mil 



Using (61) and (62), this translates into an explicit equation for Hj independent of the 
basis 

~dC^ N=1 \t,t c 



Hr 



dt 



(67) 



t = tr 



which is valid to order N = 1. Finally, we may replace C^ 4]N=1 ^ in (65) with the full 
correlation matrix C^ 4 \ The corresponding expression (67) for Hi will still be valid 
up to order TV = 1 in perturbation theory. Thus, summarizing, define the effective 
correlator 

6W(t, t c ) = C^ N=0 \t, t c )~ 1/2 C^\t, t c ) C^ N=0 \t, t c )~ 1/2 (68) 
understood as a matrix product, then the meson-meson interaction Hi satisfies 



C (4) (Mc 



-H"/(t-t c 



(69) 
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in a neighborhood of t c up to (at least) order TV = 1 in perturbation theory. 

The utility of these results in the framework of a lattice simulation lies in the 
analogy which can be drawn between (j( 4 '> N = ) an d the free correlator C^ 4 \ and 
and the full correlator C^ 4 \ The analogue of (69) may then be considered the definition 
of an effective interaction. 

B Angular Momentum Projection 

We here briefly sketch the derivation of the angular momentum-projected matrix ele- 
ments 7-6p(r,r') which appear in (47). For additional mathematical detail the reader 
is referred to ref. [6]. 

Consider the plane lattice wave \p) withp = y"(&i, k 2 ), k\ t 2 = 1 . . . L. The coordinate- 
space components are 

(f \p) = L^e** = J2 L-^Tci'JtWY^MY^fa) . (70) 

lm 

Here r is understood to be a continuous variable, with polar coordinates r and (p T , and 

Ye m (^) = -Le^, £ = 0,l...oo, (71) 

with m = if I = and m = ±1 if I > 0. Now, project out angular momentum I, m 
through 

(r{lm)\p) = d<p r Y; m (<p r ){r\p) (72) 

Jo 

= L-^J^Y^). (73) 

This defines states \r(£m)) of "good" angular momentum which have support on a 
discrete set of lattice momenta p. An operator, say K,, with continuous 0(2) symmetry 
would be diagonal in this basis 

(r(£m)\)C\s(\ri) = 6 lx 6 m ^ l \r,s). (74) 

If K, also has support on the discrete set of momenta p we may obtain its coordinate- 
space partial wave matrix elements from the expansion 

^'(v) = E E(«l(r;? 2 )X(r;p 2 )|iC|(r'; g 2 )>((r'; g 2 )|,H) (75) 

r, P 2 r',q 2 

where |(T; p 2 )) are the basis states of 0(2, Z) representations (see beginning of sec- 
tion 8). 

Now, assuming K, = Tii, we may take advantage of the 0(2, Z) invariance of Tii in 
(75) and define 

H?(r,s) = Y: E(KMI(r;p 2 ))(p 2 lM r) k 2 )((r;g 2 )|^m)}- (76) 

T p 2,q2 
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(Type) 


P 


r 


I 


f(r) 

Jim 


(0) 


Pi = P2 = 


A 1 





1/2 


(P) 


Pi > 0,p 2 = 


A 1 


4(iV- 1) 


1 






Br 


4/V - 2 


-1 


(Q) 


Pi = P2 > 


Ar 


4(iV- 1) 








B 2 


4iV - 2 


( — 1) 2 


(G) 


Pi > P2 > 


Ar 


4(iV- 1) 


cos (£7) 






A 2 
A 2 
A 2 


47V 
4/V 
4/V 


— \/2im sin(^7) 
-V2cos(^7) 
\/2im sin(^7) 



Table 4: Nonzero factors f\ m for the transformation matrix elements (77). The angle 
7 is given by tan(7) = p 2 /pi, and N £ N. 



With 7Y/ failing to be 0(2) invariant the diagonality of (74) of course also fails to 
hold. However, we may now substitute (76) for K,^(r, s) on the right hand side of (74). 
This procedure defines the rotationally invariant part of the effective interaction Tii as 
obtained from the lattice simulation. 

The transformation (r(£m)\(T; p 2 )) in (76) is easily calculated using (73) and the 
explicit expressions for |(T; p 2 )) that appear after eqn. (C.17) in appendix C of ref. [6]. 
All matrix elements have the form 

(r(£m)\(T;p 2 )) = f ^fw/ff • (77) 

For the Bessel functions Ji(x) we use the conventions of [14]. The factors are listed 
in Tab. 4. Angular momenta are parameterized as I = 4(iV — 1) = 0,4, 8... 00 or 
I = 47V — 2 = 2, 6, 10 ... 00 with N = 1, 2, 3 ... 00 in either case. A possible degeneracy 
of p 2 , which we do not encounter in the present simulation, is however accounted for 
in Tab. 4. For all combinations Y,l that do not appear in Tab. 4 the corresponding 
factors /£P are zero. 

In (76) a possible m dependence enters only through factors m 2 or im(—i) m , which 
are one in all cases. Thus m disappears completely from (76). 
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